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Abstract 

We study thermal phase transitions of color superconductivity by the lattice 
simulations of the Ginzburg-Landau (GL) effective theory. The theory is equiv- 
V/-^ ■ alent to the SUf(3) x SUc(3) Higgs model coupled to SUc(3) color gauge 

O i fields. From the eigenvalues of a 3x3 gauge-invariant diquark composite, 

a clear distinction between the 2-flavor color superconductivity (2SC) and 
the color flavor locking (CFL) phase is made in a gauge invariant manner 
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Q , The thermal transitions between the normal phase and the superconducting 

^ I phases are found to be first-order due to thermal gluons. The phase structure 

in the coupling-constant space is numerically explored and three patterns of 
phase transition, i.e., normal^2SC, normal— >CFL and normal^2SC^CFL, 
are found in the chiral limit. These results agree qualitatively with the weak- 
couphng analysis of the GL theory. 
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1 Introduction 

^ In dense and cold quark matter, color superconducting phase is expected to be realized due to diquark 

Q . pairing (see the reviews, [1|). From the analysis of the pairing interaction, the most attractive chan- 

nel is found to be in the color anti-triplet and = 0+ channel, which leads to the diquark field, 
^ ■ <Pcf = ^abc^ fjkiQ'LC'j^q'l) , as the most relevant quantity to describe color superconductivity. A three- 

^ • dimensional Ginzburg-Landau (GL) effective action written in terms of was proposed |2| to study 

D [ the thermal phase transition of the color superconductivity. Then, it was realized that the thermal gluons 

' (treated in weak-coupling approximation) turn the second-order transition in the mean-field theory to the 

first order transition |l3l|4l. Such a weak-coupling study, however, is vaUd only in extremely high density 
^ '. because the gauge coupling grows for small baryon density. Furthermore, thermal fluctuations of the 

I diquark field become also important as density decreases 1^. 

Motivated by this observation, we carry out Monte-Carlo simulations of the GL effective action 
discretized on the lattice. Although the simulation of the GL effective action is not ab initio calculation 
and is applicable only near the thermal phase boundary, it has an advantage of making non-perturbative 
calculation and simultaneously of avoiding the sign problem which is a serious obstacle in direct lattice 
QCD simulation at finite density. Since the GL action has color SUc(3) gauge invariance and the diquark 
field (Pcf is a gauge variant object, it is not a trivial task to identify the color superconducting phase in 
the lattice approach. We will see below that the eigenvalues of a gauge invariant composite ^t^x 
are useful to identify different phases, (x denotes lattice sites.) Also, analyzing the hysteresis loops of 
such gauge invariant observables, we find the first-order transitions among these phases. 

2 The Ginzburg-Landau effective action 

The 3-dimensional Ginzburg-Landau (GL) action in the chiral limit 13, which is expected to be universal 
around the thermal phase boundary, reads Sgl = / d^xC with 
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* Talk presented by M. Ohtani at the workshop on Extreme QCD (Swansea, Aug. 2-5, 2005). 
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Fig. 1: a) Phase diagram of the 3-d GL theory in the coupling constant space obtained by the mean-field theory |2|. The second 
order phase transitions are obtained, b) Phase diagram obtained by talcing into account the gluon fluctuations in a weak-coupling 
approximation. The phase boundary is shifted and the transition becomes first order. 



where we have assumed three massless quarks. Since the diquark field is anti-triplet in color space, the 
covariant derivative is defined by Di<P = {di — igT""* Af)(p. Effects of temperature and baryon chemical 
potential are incorporated in the coupling constants (m, Ai_2 and k). The action has global SUf (3) flavor 
symmetry and local SUc(3) gauge symmetry. 

Since Sgl is real, sign problem does not occur. In our actual simulation, we extend Eq.Q slightly 
to four dimension with two temporal slices and periodic boundary condition. We do so to define the 
Polyakov loop which can be used as a measure to see if the system is in the deconfinement phase. This 
extended 4-d action is expected to give qualitatively the same results as the 3-d action, although the 
parameters in the former and those in the latter (m, Ai^2 and k) are related in a non-trivial way by the 
renormalization from the temporal field-fluctuations. The GL action in Eq.Q is similar to the gauged 
Higgs model which is studied in the context of the electroweak phase transition. Note, however, that we 
have two independent coupling Ai and A2 unlike the case of the SU(2) Higgs model. 

3 Weak-coupling analysis 

The mean field approximation to the 3-d GL action leads to the second order transition between normal 
phase and the color superconducting phase |2|. Fig. la shows the phase diagram in the (Ai, A2)-plane 
in the mean field approximation. Whether 2SC or CFL is realized in color superconductor in the chiral 
limit depends on the sign of A2, and the phase boundary is located on A2 = 0. 

For weak gauge-coupling, it has been shown that the thermal fluctuation of gauge field dominates 
over the diquark fluctuation, and can be incorporated as a Gaussian fluctuation around the mean-field j 3| 
As a result, the phase transition turns into first order, which is similar to the situation in Type I metallic 
superconductor |5|. Shown in Fig. lb is a phase diagram in the (Ai, A2)-plane obtained by using the weak 
coupling approximation along the line of ref.|l3l|. Phase boundary between 2SC and CFL in Fig. la is 
shifted towards the positive A2 region by the gluon fluctuation. Furthermore, with certain parameter 
region, we observe successive transitions, normal^2SC^CFL, as temperature decreases from above. 
Since Ai = A2 is realized in the weak-coupling limit, such successive transitions may take place in high 
density QCD in the chiral limit. This possibility was first noted in |6l. 

4 Lattice simulations 

Let us discretize the 4-d version of the effective action, Eq.Q, and rescale the field and couplings as 

<P —>■ \/2a^/o, Aj — > Xijc?, rn^ ^ (1 — 2Ai — 2A2 — ^R)/aa?, k —>■ R/a (a being the lattice constant). 
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Fig. 2: Eigenvalues of 'I't'^x as a function of update steps. (Ai, A2) — (0.03, 0.01), R = 0.34, /3 — 5.1 for the left panel 
a) and (Ai, A2) = (0.03, 0.03), R — 0.43, /3 = 5.1 for b). The isolated large eigenvalue associated with two degenerate small 
eigenvalues a) corresponds to 2SC phase while three degenerate eigenvalues after thermalization b) indicate CFL phase. 



Then, we reach the following lattice action 



S = Y,\ Tr(4^x) + (Ai + A2) (tV«.^,) - 1 j + A2 ( Tr(^t 



(2) 



where Sg is the standard plaquette action with /? = To ensure that the system is in the deconfined 
phase, we keep (3 around 5.1 according to the critical value obtained in the pure Yang-Mills case with 
Nt = 2. Taking several sets of the quartic couplings Ai 2, we make simulations by scanning the param- 
eter K. In this article, we report the results obtained in the simulation with a lattice volume of 2 x 12^. 
We note that the 2x8'^ lattice gives qualitatively the same results. We use the pseudo heat-bath method 
to update the gauge-hnk variables and generalize the efficient update-algorithm of SU(2) Higgs-field ITJ 
to our case. 

The color superconducting phase and the normal phase are distinguished by looking at the gauge 
invariant operator Tr ^l:<Px- Although the field fluctuations give non-vanishing expectation value of 
this operator even in the normal phase, we can observe discontinuous changes of the expectation value 
as K changes. In addition to this, we measure other gauge-invariant operators such as Tr|D^(^p, 5g and 
the Polyakov loop to confirm that the discontinuous change takes place for all these operators at the same 
K. (Note that this discontinuous change of the Polyakov loop does not mean deconfinement transition.) 
This constitutes an evidence of the first order phase transition. 

Moreover, we can differentiate different color superconducting states (2SC, CFL or something 
else) by diagonalizing the 3x3 matrix ^^t^x in the flavor space and extracting its eigenvalues. As 
shown in Fig.2a, we observe one large eigenvalue with two degenerate small eigenvalues in a particular 
set of (Ai, A2). In this case, we regard the state as 2SC in which the quark pairing takes place in one 
specific color-flavor channel. We also observe three degenerate eigenvalues as shown in Fig.2b with 
another parameter set, which corresponds to the CFL phase where all channels are involved equally to 
the pairing. Here we emphasize that this is a gauge invariant way of identifying different phases in color 
superconductivity. As we scan the parameter R, we find rapid changes of the eigenvalues of ^t^x 
at critical value R^. When we simulate around this point by increasing and decreasing R, we observe 
hysteresis loops in the eigenvalues as well as in the Polyakov loop, as shown in Fig.3. Such observation 
also bears evidence of the first order transition. 

Near the critical point of the first order transition, it is always difficult to find a global minimum 
of the free energy among nearly degenerate local minima. To find a true minimum, we divide the lattice 
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Fig. 3: The hysteresis loop of J2x Tr<P!c<Px obtained by increase or decrease of k. (Ai, A2) = (0.01, 0.005) and /3 = 5.1. 

volume into domains so that we can put different phases obtained by the hysteresis loop into different 
domains. Starting from such a mixed configuration, a state with the highest pressure pushes away the 
other metastable states after many update steps. We call this as 'boundary-shift test'. In this way, we can 
determine the most favorable state for each value of R, and can estimate Re- 

In Fig.4, we show the phase diagram in the (Ai, A2)-plane obtained by our lattice simulations. The 
points marked by the symbols (star, cross, square, and plus) correspond to the places of actual numerical 
simulations. We find three types of thermal transitions depending on the coupling parameters: for large 
values of A2, transition from the normal phase to CFL takes place as star and cross symbols show. For 
small or negative values of A2, transition from the normal phase to 2SC takes place as plus symbols show. 
In between the two cases, there is a parameter region where successive transition, normal— >^2SC—>^CFL, 
takes place as square symbols show. 

The global feature of the phase diagram in Fig.4 is qualitatively similar to that of the weak- 
coupling analysis shown in Fig. lb, although the coupling constants (e.g. Aj vs. Aj) are related in a 
nontrivial manner. Note that, near the phase boundary between 2SC and CFL, we find a region (shown 
by the cross symbols) where metastable 2SC survives, even when the global minimum of the free energy 
transfers directly from the normal phase to CFL. This metastable 2SC state, which is a local minimum 
of the free energy, can be seen in the hysteresis loop but is pushed away by the 'boundary-shift test' . The 
jump of the operators at the transition points becomes smaller if we take larger Aj. It is considered that 
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Fig. 4: Phase diagram in the coupling constant space ob- 
tained by the lattice simulation. Symbols show the pa- 
rameters at which the actual simulations are carried out. 
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Fig. 5: Contour plot of the free energy at the transition 
point from normal phase to CFL in the weak coupling 
analysis. Scales are in an arbitrary unit. 



5 



thermal fluctuations of the diquark field weaken the first order transition. 

Fig.5 shows the contour plot of the free energy at a critical point of the transition from the nor- 
mal phase to the CFL phase obtained in the weak-coupling approach. Setting the diquark field as 
diag(Ai, Ai, A3), the gluonic fluctuations are calculated as a function of Ai and A3. The dark re- 
gions have smaller free energy than the light regions. The normal, 2SC and CFL phases are identified 
as Ai = A3 = 0, (Ai = 0, A3 / 0) and Ai = A3 7^ 0, respectively. This figure indicates a subtle 
interplay among three local minima as a function of temperature and coupling constants. It also serves 
as an possible explanation of the phase structure shown in Fig.4. 

5 Discussion 

We have studied thermal transitions of color superconductivity on the basis of the lattice simulations 
of the Ginzburg-Landau effective action. The relevant degrees of freedom are the diquark field (p and 
the gluon field. We have identified different phases of color superconductivity by measuring the gauge 
invariant flavor-matrix J2x '^x'^x and extracting its eigenvalues. 2SC and CFL phases as well as normal 
phase are clearly distinguished in a gauge invariant way. 

In our simulations, we observed the hysteresis loop of various gauge invariant operators as a 
function of the parameter R, which designates the first-order thermal transition. The global minimum of 
the free energy for each k is determined by the boundary-shift method in which we carry out updates by 
starting from a mixed configuration obtained in the hysteresis loop. With these procedures, we obtained 
the phase diagram in the coupling constant space and compared the result with that of the weak coupling 
calculation. Three types of thermal transitions are observed depending on the couplings: normal^2SC, 
normal^CFL and a novel transition, normal^2SC— >CFL. 

To find a realistic phase diagram of color superconductivity in the temperature (T) and chemical 
potential (//) plane, it is necessary to make a connection of our coupling constants to T and fi, which 
requires a matching of the GL effective theory to the microscopic QCD. The quaik masses and color 
and charge neutralities become important at low baryon density and will lead to more complex phase 
structure than that shown in this report. Further studies are needed to attack these problems. 
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